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ABSTRACT 

In this contribution, we test our previously published one-dimensional PDR model for 
deriving total hydrogen volume densities from HI column density measurements in ex- 
tragalactic regions by applying it to the Taurus molecular cloud, where its predictions 
can be compared to available data. Also, we make the first direct detailed comparison 
of our model to CO(l-O) and far-infrared emission. 

Using an incident UV flux Go of 4.25 (x = 5) throughout the main body of 
the cloud, we derive total hydrogen volume densities of ~ 430 cm"'^, consistent with 
the extensive literature available on Taurus. The distribution of the volume densities 
shows a log-normal shape with a hint of a power-law shape on the high density end. 
We convert our volume densities to H2 column densities assuming a cloud depth of 5 
parsec and compare these column densities to observed CO emission. We find a slope 
equivalent to a CO conversion factor relation that is on the low end of reported values 
for this factor in the literature (0.9 x 10^" cm^^(K km s~^)~^), although this value is 
directly proportional to our assumed value of Go as well as the cloud depth. We seem 
to under-predict the total hydrogen gas as compared to 100 /im dust emission, which 
we speculate may be caused by a higher actual Go incident on the Taurus cloud than 
is generally assumed. 
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1 INTRODUCTION 

Molecular hydrogen is widely considered the key ingredi- 
ent in star formation. However, on extragalactic scales it is 
only observed directly at higher densities or excitation condi- 
tions. Photodissociation regions (where molecular hydrogen 
is destroyed) are key to learning more about the molecu- 
lar gas, because observable amounts of atomic hydrogen are 
produced in these regions. 

The available literature on photodissociation regions 
(PDRs) is extensive. Early work includes HoUenbach et al 



( 1971 1 and Jura ( 1974 1 on the formation and destruction of 
molecular hydrogen, [Glassgold fc Langer| ( |1975[ ) on the ap- 
pearance of C"'' and CO, and calculations on more atomic 
and molecular species by e.g. Black & Dalgarno ( 1977 1 and 



van Dishoeck & Black ( 1986 1. Recent years have seen the ar- 



rival of various publicly available computer models that aim 
to solve the detailed structure of PDRs for many atomic and 
molecular species, such as CLOUDY ([Ferland et al.|p98l ), 
the Me udon PDR code ([Le Petit et al.|2006[ ), the K OSMA- 
T code pBllig et al.||2006| ) and the 3D-PDR code ( [Bisbas 
et al.|2012 I. A variety of such codes were compared and the 
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results published by RoUig et al. (20071. For the focus of this 
work, we will limit our attention to atomic and molecular 
hydrogen (as well as, to a certain extent, carbon monoxide) 
in low density molecular clouds, since we aim to explore and 
provide an independent probe of total hydrogen volume den- 
sities as opposed to a detailed model of the internal structure 
that includes predicted line intensities of various species. 

In lower-density, lower-radiation regimes, it is possi- 
ble to derive total hydrogen densities using the relatively 
straightforward process of photodissociation of molecular 
hydrogen, which results in detectable atomic hydrogen emis- 
sion, potentially on galactic scales. A detailed study of a 
one-dimensional approach to such photodissociation regions 
was carried out by |Hollenbach et al.| ( |199T| ), including calcu- 
lations connecting gas densities with the incident UV flux, 
CO and 100 /xm dust emission. Their results have been used 
to model the interstellar medium of external galaxies as well 
as the Milky Way. 

|Li et al.| ( |2002| ) studied PDRs in what they called the in- 
termediate regime, or n « 500—5000 cm~'^ and Go ~ 10—60, 
in the L1204 and S140 regions. While they successfully pre- 
dicted several atomic and molecular emission lines, they 
observed enhanced H2 emission that could possibly be ex- 
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plained by enhanced formation rates or weak shocks related 
to turbulent decay. 

The notion that a large fraction of extragalactic atomic 
hydrogen in galactic disks is produced in PDRs has been 



around for quite some time ( 


Allen et al.||1986| |1997| |Smith| 


et al.|2000 Heiner et al.|2011 


I. Using PDR-produced atomic 



hydrogen, it becomes possible to calculate total hydrogen 
volume densities at selected positions across galactic disks 
using a simple approximate analytical one-dimensional PDR 
model (see also e.g. Higdon et al.|20lT l. The resulting den- 



sities are consistent with values expected for the diffuse in- 
terstellar medium. Corresponding molecular emission (CO) 
is found near higher (several hundred cm"'') total hydrogen 
volume densities, although it should be noted that under- 
excited CO could still hide in relatively dense clouds. The 
densities obtained in this manner are also consistent with 
the (volumetric) Schmidt Law, linking the local star forma- 



tion rate to the total hydrogen volume density ( Heiner et al 



2010), at least in the power-law exponent. 

However, so far this approach has not been compared 
directly to CO-derived molecular hydrogen gas contents in 
a detailed fashion. In order to do this, we turn to our own 
Galaxy to apply the one-dimensional PDR model to a molec- 
ular cloud that has densities in the range where this model 
can be expected to apply, namely below 10^ cm~^. 

We selected the well-known and heavily-studied Taurus 
molecular cloud, part of the Perseus- Taurus- Auriga complex 
for this study. It has molecular gas densities of a few hundred 
cm~^ with denser clumps (for a recent comprehensive study, 
see e.g. Pineda et al. 20101 and is relatively quiescent in 
terms of star formation. 

The atomic hydrogen in the Taurus area is considered 



to be connected to the molecular gas ( Li & Goldsmith 2003 1. 
Several surveys of Taurus have been completed in the CO 1-0 
line ( [Ungerechts fc Thaddcus 1987: Narayanan et al.||2008| ) 
as well as a number of studies focusing on specific regions. 



According to Ungerechts & Thaddeus (19871, the Tau- 
rus CO mass is of the order of 10 M0, depending on which 
parts of the complex are included (see also [Pineda et al.| 
20101. Furthermore, clouds in the Taurus complex were 



found to have CO-derived masses approximately equal to 
their virial masses ( |Heyer et al.|1987| . [Ab"ergel et al.| ( [19951 ) 
compared the Taurus CO emission with IRAS far-infrared 
emission maps to conclude that the CO emission should 
trace at least 70% of the matter indicated by the 100/im 



emission, using the Ungerechts & Thaddeus ( 1987 1 maps 



This makes Taurus a good candidate to compare the results 
of the simple one-dimensional PDR model to known molecu- 
lar data, which will support the case for applying this model 
in the extragalactic case as well. 

In the case of Taurus, then, we have a molecular cloud 
that is bathed in a photodissociating flux coming from the 
outside. We may assume that the molecular cloud is sur- 
rounded by an envelope of atomic hydrogen, that is either in 
the process of becoming molecular, or originates from pho- 
todissociated molecular hydrogen. Either way, this atomic 
hydrogen is part of the balance of photodissociation. Our 
study provides a starting point to further investigation of 
the atomic-to-molecular transition around molecular clouds 
in the context of the kinematics of the g prompted by, 
e.g., the work presented b y[Ballesteros-Paredes et ar] ( [1999| 
and Nguyen Luong et al. (2011 1. 



The application of the simple one-dimensional PDR 
model includes an assumption about the details of the con- 
version of HI to H2 and back, namely that (at least in the 
extragalactic case) the full integrated HI column can be as- 
sumed to be part of the process, minus a local background 
value. For an edge-on large-scale PDR, this assumption is 
likely to overestimate the HI column and therefore underes- 
timate the related total hydrogen volume density. The rea- 
son for this is that the line of sight from the observer through 
the galaxy probes edge-on HI shells more or less parallel to 
the shell's surface, which would yield a higher column den- 
sity than the column density seen at a right angle to the 
surface (assuming a shell that is thinner than it is wide). 
The true three-dimensional structure of the cloud remains 
unknown, which is especially true for the extragalactic case. 

In our own Galaxy, Allen et al. (2004) used a limited 



velocity range dictated by the local molecular emission (for 
Maddalena's Cloud), an approach that we also adopt here. 
Under these assumptions, we can trivially calculate the frac- 
tion of the total HI column density that is included in the 
chosen velocity range. However, this in itself only hints at 
the fraction of HI that is part of the transition from atomic 
to molecular gas (and back). 



Also, while Allen et al. ( 2004 1 applied the simple PDR 



model to Maddalena's Cloud in our own Galaxy, a detailed, 
spatially resolved application fell outside the scope of their 
work. With the application of this method to Taurus, then, 
we expand on their result. Consequently, this allows us to 
study the distribution of total hydrogen volume densities in 
the Taurus region. This distribution holds important clues 
to the dynamical state of the cloud. 

This paper is structured as follows: First, we briefly re- 
view the theory behind the one-dimensional PDR model and 
how it can be applied to the case of Taurus. We then describe 
the archival data we used to explore the Taurus region. We 
discuss obtaining appropriate atomic hydrogen column den- 
sities, after which we present total hydrogen volume density 
maps of the main Taurus area. We then proceed to com- 
pare our results to CO emission data, as well as far-infrared 
data. We check our results for consistency with the |Allen| 
et al. (2004) calculations and discuss the values of the UV 



field suggested by our study. Finally we end with a summary 
of our conclusions. 



2 A SIMPLE ONE-DIMENSIONAL PDR 

MODEL APPLIED TO THE MOLECULAR 
CLOUD-ATOMIC ENVELOPE INTERFACE 

We use a straightforward one-dimensional plane-parallel 
PDR model that connects atomic hydrogen, molecular hy- 
drogen (through the total hydrogen volume density), the 
incident UV flux and the local dust-to-gas ratio. For a more 
thorough description we refer to, e.g. 



Heiner et al. ( 2008 ) 



The flnal, approximate formula describing the photodissoci- 
ation balance, which we employ here, is: 



n = 106 Go 



exp 



7.8 X 1020 ) 



(1) 



where n is the total hydrogen volume density in cm"'' aver- 
aged along the line of sight through the cloud. Go is the inci- 
dent UV flux and S/Sq is the dust-to-gas ratio scaled to the 
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solar neighborhood value and A''hi is the (PDR-produced) 
HI column density in the direction of the incident UV flux. 
The exponent in 5 /So is slightly difi'erent from the one that 
was used previously because it includes the improvements 



suggested by Heaton ( 2009 1 



The atomic hydrogen column density predicted by this 
approximate analytical model is scaled to fit the column den- 
sities produced by numerical PDR calculations presented by 



Kaufman et al. (19991, as described in Allen et al. (20041. 
It assumes that the formation and destruction of molecu- 
lar hydrogen are balanced, at least when integrated through 
the gas cloud (Sternberg 1988). It provides an estimate of 
the total hydrogen volume density at the location where a 
certain HI column density is measured, neglecting the finer 
details of the formation and destruction of molecular hy- 
drogen and other elements in order to focus on the direct, 
line-of-sight-averaged connection between atomic hydrogen 
and molecular hydrogen. 

In the current context. Go = 0.85x ( [HoUenbach fc Tie-| 
lens 1 1999 1, where Go is the incident radiation in units of the 
Habi ng flux |Habing|[l968t and X is in units of the Draine 
fleld ( |Draine| 1978[ ). These definitions of the radiation field 
and their relation are valid under the assumption of a cloud 
approximated by a semi-infinite slab, illuminated isotropi- 
cally, with radiation therefore entering through 2ti steradi- 
ans. In the case of Taurus, it is not clear where the incident 
flux originates, and we approximate the ambient radiation 
fleld by adopting a constant value for it following [Goldsmith] 
et al. (20101. Equation ([TJ requires a single value for the in- 
cident flux, namely the one at the edge of the cloud. In our 
model, it is not necessary to know whether the radiation 
originates from in front of, or from behind the cloud. The 
value of n, however, is an average and includes all the hy- 
drogen gas along the line of sight through the cloud. This 
includes the atomic hydrogen at the edge of the cloud as 
well as the molecular hydrogen further inwards. This aver- 
age value may include clumped gas along the line of sight as 
well. 

This model is thought to be most useful in th e den- 
sity range 10 cm"^ < n < lO'^ cm"^ and Go < 100 ( [Allen 
et al.|2004 l. The underlying assumption is that (a large frac- 



tion of the) measured HI column density is indeed produced 
by photodissociation due to the nearby source of UV flux. 
Trivially, the total hydrogen column density N = n x dcioud 
where ddoud is the assumed cloud depth. Throughout this 
paper, n or N without subscript will refer to total hydrogen 
quantities, that is A'^ = Ahi -I- 2Ah2 ■ 

In the extragalactic case, one needs to identify patches 
of HI in column density maps that are suspected to be as- 
sociated with central OB star associations as seen on far- 
UV images. The linear dimension of the PDR model is then 
taken to stretch on the plane of the sky from the central star 
cluster through the HI patch. From the projected angular 
separation, a linear distance and a corresponding incident 
flux on the location of the HI emission are then calculated. 
It is important to select a PDR that is seen approximately 
edge-on in order to be able to calculate the incident flux, 
although we will suffer from projection effects. While the 
HI column density is not measured in the direction from the 
star cluster through the HI patch, it is nevertheless assumed 
that the measured column density is an appropriate probe 
of the required column density, as we explained in the In- 



Radial velocity (km/s) 

Figure 2. Sample profile from the GALFA HI datacube of Tau- 
rus. This particular profile shows a clear absorption feature as 
shown by three gaussians used to fit the line. Thick solid lines; 
HI profile and CO profile (multiplied by a factor of 5 for clarity). 
Dash-dotted lines: gaussian components of the fit. Dashed lines: 
fitted profile and residual. Gray area: section of the profile used 
to infer the local HI column density. 



troduction. Specifically, we refer the reader to |Heiner et al.| 
(20081, their Figure 4c, showing an idealized representation 



of the expected morphology. Formally, a spherical molecular 
cloud is assumed, or a tubular filament in the plane of the 
sky. A more thorough study of the projection effects on the 
measured HI column density in the extragalactic case falls 
outside the scope of this paper, and will be investigated in 
more detail in a future contribution. 

The applicability of the method is limited by the angu- 
lar resolution of the images, meaning that only large-scale 
PDRs can be used (about 50-500 parsec separation). To ap- 
ply this same method to PDRs in our own Galaxy, we take a 
slightly different approach. Instead of considering an edge- 
on PDR, we look at the PDR face-on. Here, we consider 
a uniform radiation field impinging on a known molecular 
cloud along the line of sight (with the radiation coming ei- 
ther from the direction of the observer or from behind the 
cloud), and attempt to isolate the photodissociated atomic 
hydrogen in front and behind the molecular cloud. The lin- 
ear dimension of the PDR model is therefore perpendicular 
to the plane of the sky, and we can consider each pixel on 
the observed map individually, as long as the approximation 
of the incident radiation going along the line of sight holds. 
In principle, we should now measure the HI column den- 
sity in the proper direction, although the three dimensional 
structure of the cloud remains unknown. 

In this paper, we will compare our predicted column 
density directly to the observed (velocity-integrated) 

Wco CO(l-O) line emission. This CO emission is rou- 
tinely converted to H2 column densities using a constant 
conversion factor, otherwise known as Aco- In principle 
this quantity is defined only for an ensemble of molec- 
ular clouds on a galactic scale, but this conversion has 



also been applied on sub-kpc scales (Bigiel et al 



or specific regions of our Galaxy (Abdo et al. 20101. The 



2008 1 



value of this so-called X-factor is thought to range be- 
tween 0.9-4 X 10^'' cm"2 (K km s"^)"\ but a Galac- 
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CO 4.BB km/s; HI 4.73 km/s 



CO 6.18 km/i; HI 6.25 km/s 




CO 6.83 km/Si HI 6.99 km/s 



CO B.13 km/s; HI 7.73 km/s 




4" 40'^' 4" 20' 

R.A. (2000.0) 



Figure 1. HI and CO emission near 5, 6, 7 and 8 km s . CO contours are main beam brightness temperatures T^i, of 1, 3.5 and 6 K. 



tic standard value in the same units is usually adopted for 
any galaxy in the absence of more direct measurements of 



this value. For example, Kennicutt ( 1998 1 used a factor of 
2.8 X 10'^'^ cm~'^ (K km s~ )~ for a sample of 61 galaxies 
to derive H2 column densities. The lower values have come 
from e.g. Digel et al.| ( |1996 1 and the higher values from e.g. 



Draine et al. ( 2007 1 



For the purpose of our comparisons we will use the CO 



X-factor value adopted by Dame et al. (2001 1, namely 1.8 x 
lO'^" cm~'^ (K km s~^)~'^, since we use their CO data in this 
work. It should be noted that we use this value for reference 
purposes only, as it is helpful to plot this value as the slope 
in our figures comparing to Wco- At our parsec and 
sub-cloud scale a constant conversion factor should not be 
expected. 



3 DATA AND PARAMETERS 

For the one-dimensional PDR model described in Equation 
|l]) we need the following parameters: the incident UV flux 
Go, HI column densities, and the dust-to-gas ratio (scaled 
to the solar neighborhood) 5 /So- In addition we will be com- 
paring our predicted densities to CO(l-O) and far-infrared 
(100 fim) derived H2 column densities. 



incident x of 5 (or Go = 4.25). They consider this value on 
the high end of the range based on a census of nearby B- 
stars, but they require at least a value this high to explain 
their observations. In our model, n oc Go, which means that 
underestimating Gq also means underestimating the total 
hydrogen volume density. 

We assume a solar dust-to-gas ratio (5/ So — 1). Fixing 
both the incident UV flux Go and the dust-to-gas ratio, ap- 
propriate for a relatively small region like the Taurus molec- 
ular cloud (as opposed to large-scale PDRs up to several 
hundred parsec across as investigated in the extragalactic 
case), means that the predicted total hydrogen density will 
vary with the HI column density only. 

We use the 21-cm data from the GALFA-HI survey 



which has an effective beam size of 3.9 x 4.1 arcmin (Peek 
|et al.|[2bll[ ). We retrieved the data with a velocity separa- 
tion of 736 m s~^, ranging from a radial velocity of -10 km 
s~^ to 20 km s~^, which is more than enough to include the 
Taurus molecular cloud at a radial velocity between 5 and 
6 km s~^. 



Ballesteros-Paredes et al. ( 2009 1 modeled the main part 



of the molecular cloud as a 32 x 5 pc fllament at a distance 
of about 145 pc, based on the distance measurements of 



Following Goldsmith et al. (20101, we adopt a uniform 



Loinard et al. (20071 and Torres et al. (20071. We wiU adopt 
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HINSA 5,26 km/s; CO 6.16 km/5 

1 7.0 




4i'4Q". 41.20'" 4'' 40'" 4'' 20"' 

R.A. (2000.0) R.A. (2000.0) 



Figure 3. Map of the second derivative of the HI profiles indica- 
tive of the strength of HINSA features (integrated over a range 
of 5 to 7 km s"-*^). 



Figure 4. Contours of HINSA strength near 6 km s ^ on top 
of CO emission in grayscale. HINSA strength contours are 4, 12 
and 20 K km-2 s^ 



this model here. The actual distance does not influence our 
results. 

We use the CO data of Taurus presented by |Dame et al.| 



( 2001 1 for comparison. We reprojected it to the HI map first, 



but this results in an over-sampled CO map. We therefore 
smoothed the HI data to the CO resolution of 8.4 arcmin. 
Then, we regridded the data to a common grid; from 1 ar- 
cmin to 8 arcmin for the HI data by decimation, and from 
7.5 arcmin to 8 arcmin by reprojection in the case of the CO 
data. This allows us to compare the data on a pixel-by-pixel 
basis. 



Wolfire et al. ( 1993 1 investigated the expected CO line 
emission from GMCs and found that, for low UV fluxes 
(1 < Go < 100), the CO luminosity varies by only 10%. 
In order for the CO line profiles to appear smooth and cen- 
trally peaked, a large number of clumps need to be included 
in the (observing) beam. Since the Taurus CO line profiles 
indeed show this smooth shape, we believe assuming a uni- 
form Go is justified. 

We also used the 100 /im far-infrared dust map 



(temperature-corrected) produced by Schlegel et al. ( 1998 1 



by extracting the general Taurus area from their all-sky map. 
The angular resolution is 6.1 arcmin, which we smoothed 
and regridded to 8 arcmin to compare the results. 



4 OVERVIEW OF THE TAURUS REGION 

In Figure [l] we show a number of CO/HI overlays that il- 
lustrate how the atomic and molecular gas in the Taurus 
area are intimately connected. Qualitatively, there is an anti- 
correlation between the CO and HI emission: the HI emis- 
sion is clearly depressed in the areas where CO occurs. We 
aimed to locate parts of the Taurus cloud where our one- 
dimensional PDR model is likely to be applicable. In other 
words, we needed to find the parts on the map that seem 
to feature CO emission surrounded at the front and at the 



back by HI emission, along our line of sight. We then as- 
sumed that this HI emission is related to the balance of 
photodissociation described in the simple PDR model. 

We show a fairly typical HI profile with the correspond- 
ing CO profile (artificially magnified by a factor of 5 for clar- 
ity) in an area dominated by CO emission in Figure [2] At 
first glance, this profile seems to show either two overlapping 
HI emission peaks, or perhaps a single peak with an absorp- 
tion dip at the center. Both scenarios may point to the con- 
figuration we aimed to find: atomic hydrogen surrounding 
carbon monoxide, or atomic hydrogen tracing dense gas by 
displaying a (self-) absorption feature respectively. A good 
fit of this particular profile can be obtained using three 
gaussians (shown in the figure): a broad peak, small nar- 
row peak and one gaussian representing absorption. Poorer 
fits are obtained using two or three positive-only gaussians 
(not shown) . While many pixels in the datacube of the Tau- 
rus area can be fitted through this scheme, this in itself does 
not necessarily imply that absorption indeed occurs. 

|Li fc Goldsmith] ( |2003[ ) conducted a survey of dark 
clouds in the Taurus/Perseus region in search of narrow 
absorption features in HI and OH emission. If the HI ab- 
sorption feature is narrower than the CO linewidth, they 
call it HI narrow self-absorption (HINSA). They found that 
HINSA features correlate well with molecular tracers. For 
the Taurus area, they found HINSA column densities of 
about several times 10^* cm~^, where it should be noted that 
measuring the HINSA feature does not necessarily mean 
fully filling the 'dip' in the line profile. 

Indeed, the (negative) absorption gaussian profile in our 
Figure [2] corresponds to a column density about 30 times 
higher than 1 x 10^* cm~^, which would lead us to over- 
estimate not just the absorption, but also the HI emission 
grossly if we were to add both the positive broad gaussian 
as well as the absorption gaussian to infer a total HI col- 
umn density. Therefore, we decided against fitting the HI 
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Figure 5. Fraction of the HI column density contained in our 
chosen velocity range of 4 — 8.5kms~^ relative to the total HI col- 
umn density over the range of -10 to 20 km s"'^ (both background- 
subtracted) . 



column fractional error 
overestimation of n 



Distance from fiiament center (pc) 

Figure 6. Estimated errors when the HI column is not taken ex- 
actly through the center of a filament, assuming a diameter of 5 
parsec. Dashed line: fractional error in the HI column, at the fil- 
ament center and slowly increasing going outward. Dashed-dotted 
line: (fractional) overestimation of the calculated total hydrogen 
volume density. It slowly increases up to 50% at about 2 parsec 
away from the center and then starts to increase rapidly. 



line profiles in tlie Taurus region as a way to identify areas 
wliere the simple PDR model can be applied. 

Considering the comparatively low HINSA column den- 
sities found by Li & Goldsmith ( 2003 ) , and in a more physi- 



cally refined way by |Krco et aL[ ( |2008| , relative to the overall 
HI column density, we finally elected to ignore the HINSA 
features for the purpose of calculating the HI column density 
related to the PDR. Instead, we opted to use the presence 
of HINSA only as an indication that our simple PDR model 
applies, independent of the presence of any CO emission. 

As mentioned in the previous Section, we assume that 
the main Taurus molecular cloud can be pictured as a 32-by- 
5 pc filament. We can use this morphological information, 
specifically the filament depth, to estimate a total hydrogen 
column density from the volume densities obtained by our 
method. Choosing the appropriate line of sight is important, 
since due to the approximative nature of the model, it is 
possible for it to arrive at total hydrogen column densities 
that are less than the HI column density. It is therefore 
essential to make sure the model is applied to the proper 
regions. In the extragalactic case this is not a problem, since 
the columns are sufficiently deep not to get this close to 
the regime where the one-dimensional PDR approximation 
breaks down. 



5 ATOMIC HYDROGEN COLUMN DENSITY 

As mentioned in the previous section, even though we do 
not account for HINSA column densities in our HI calcu- 
lations, we do use the presence of HINSA features to trace 
the pixels where we expect the simple PDR model to apply. 
It should be noted however that this does not mean that 
the model does not apply elsewhere in the cloud, so we do 
not necessarily trace all parts of the cloud where our model 
is applicable. To find HINSA, we use the method described 



in Krco et al. (20081, where they show that HINSA is most 



easily seen in the second derivative of the HI line profile with 
respect to the velocity. The reason for this is that HINSA 
absorption features are typically Gaussian-like. Derivative 
functions of a Gaussian show a progressively stronger de- 
pendence on line width, making the narrow HINSA feature 
an easy-to-identify dominant feature in the second deriva- 
tive as the broader overall profile is smoothed out. We can 
therefore construct a map that is indicative of the strength 
of the HINSA feature across the Taurus region, as shown in 
Figure |3] We integrated the second derivative velocity chan- 
nels (in units of K km"'^ s'^) over the range of 5 through 7 
km s"^ Then, we chose to select only those pixels with an 
integrated value of 4 K km~^ s or higher. We will refer to 
this integrated value as the HINSA strength. The resulting 
map shows the main Taurus area as well as features in other 
parts of the region. Figure |4] shows how the HINSA strength 
traces part of the Taurus CO emission. 

The calculation of the atomic hydrogen column density 
(ignoring HINSA) requires additional assumptions. Given 
the very large width of the HI line (~ 20 km s~^; Figure [2]), 
it is natural to conclude that it originates from gas extending 
over much larger distances along the line of sight than the 5- 
pc assumed thickness for the cloud. Thus, in order to select 
only the atomic gas directly participating in the atomic to 
molecular transition into the cloud, we integrate only over 
the same velocity range that is spanned by the molecular 
gas at the position of the cloud, around 6 km s~^, from 4.0 
km s~^ to 8.5 km s~^ (gray area in Figure|2]). For our back- 
ground subtraction we simply took the minimum value in 
each profile over the range of -10 to 20 km s~^, and we con- 
verted the observed velocity-integrated brightness temper- 
ature Tb to A^'hi in the usual optically thick approximation 
by multiplying by 1.823 x 10^* cm"^ (K km s^^)"\ Con- 
sidering that Taurus is located at the Galactic anti-center, 
we do not expect any noticeable effects of galactic rotation 
on the radial velocities. A similar strategy to isolate the 
atomic hydrogen connected to a molecular cloud was used 
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Figure 7. Total hydrogen volume densities in the main Taurus area, full (left) and decimated (right) version (used to compare with 
CO emission). 



by e.g . Allen et al. ( 2004 1 and for Taurus itself by Hurwitz 
( 19941, although the latter authors used a 4 km s~^ window 
centered on a 4 km s~^ radial velocity. Their HI column den- 
sities, only derived at selected positions, range from 1.6 to 
2.7 xlO^'' cm~^. Despite the different central velocity, this 
is comparable to the low end of our values, which also in- 
clude the higher HI column densities further away from the 
main Taurus region. Figure [5] shows the fraction of the total 
HI column density (background-subtracted) we recover in 
our chosen velocity range as compared to integrating over 
the full velocity range. Typically, this fraction is between 0.3 
and 0.4. 



We have also estimated the errors we make in the col- 
umn and volume densities when assuming to be looking 
straight through a tubular filament. Our model assumes that 
the HI column is measured along the incident flux and, since 
we assume a uniform field impinging onto the surface of a 
filamentary (cylindrical) cloud, we should always measure 
the HI column straight through the center of the filament. 
In practice, we have a flat image and we do not know the 
exact fllament center, so an error will be made when mea- 
suring the HI column directly from the map. Figure [6] shows 
the error made by taking the HI column off-axis, purely as 
a geometrical effect due to underestimating the path length 
through the filament. Obviously, the error close to the edge 
of the cloud is the largest, but the error in assumed HI col- 
umn density is relatively modest for most of the distance 
between the center and the edge of the cloud. At 2 parsecs 
from the center (80% of the way) we would underestimate 
the HI column by about 50% and overestimate the total hy- 
drogen volume density by about a factor two. Considering 
that the HINSA features seem to occur well inside the cloud 
(on the plane of the sky) and not too close to the edges of 
the main Taurus area, we conclude that at least in the main 
part of the molecular cloud these errors do not dominate our 
results. 



6 TOTAL HYDROGEN VOLUME DENSITY 
MAPS 

We calculated a total hydrogen volume density map at the 
native HI resolution and at a reduced resolution of 8 arcmin 
used to match the CO data. The full resolution map and the 
decimated map are shown in Figure [7] We stress that this 
calculation could be performed over a larger region. How- 
ever, our model only works when we measure HI column 
densities close to the axis of the cloud filament - for lines of 
sight approximately perpendicular to the cloud surface. We 
assumed that cloud filaments are traced approximately by 
the level of HINSA strength. 

The distribution of the calculated total hydrogen vol- 
ume densities is shown in the histograms in Figure [8] In the 
left panel, the total hydrogen volume densities are shown 
for all pixels having a HINSA strength of at least 4. We 
find an average volume density of 436 cm~'^ and a standard 
deviation a — 129 cm""^. If we look at specific regions of 
the cloud with strong CO emission (next section) we obtain 
a more limited range of values (right panel), with densi- 
ties around 500 cm~^. A total hydrogen volume density of 
about 430 cm~"^ is consistent with the 405 cm"'' assumed 
by Ballesteros-Paredes et al. (20091, although theirs is an 



average for the whole cloud and does not include atomic 
hydrogen. 

The shape of the density distributions is expected to 
be log-normal for isothermal fiows, with a power-law tail for 



deviations due to gravitational collapse ( Vazquez-Semadeni 
1994; Pas sot fc Vazquez-Semadeni|1998l 



Vazquez-Semadeni 



jetaL2008l. More specifically, the shape of the density dis- 



tribution should evolve based on the dynamics of the cloud, 
and the exact value of the power-law slope depends on the 



projection of the cloud against the sky as well (Kainulainen 



et al. I [20091 |Ballesteros-Paredes et"ar]|20TT| [Kritsuk et al 
2011 1. The left panel of Figure |8| shows a simple log- normal 
fit to the distribution. There is a hint of a power-law tail at 
the high-density end, whose existence would indicate that 
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Figure 8. Left panel: Distribution of logarithmic total hydrogen volume densities using a HINSA strength threshold of 4., with a simple 
fit of a gaussian distribution. Note that this histogram is taken from the full-resolution density map. Right panel: Same, but only for 
regions 2, 3, 4, and 5 defined in Section [7] 



the high density gas of the cloud is dominated by gravity. A 
more sophisticated determination of the HI column density 
as well as a more detailed knowledge of the radiation field 
Go would improve the analysis of the shape of the density 
distribution. 

Again assuming an ellipsoidal shape for this region of 
Taurus, we arrive at a total hydrogen mass of ~ 1 x 10*Mq, 
consistent with the CO mass in |Ungerechts fc Thaddeus] 
(19871). 



7 COMPARISON TO CO EMISSION 

Assuming a cloud depth gives us the total hydrogen column 
density from the total hydrogen volume density and allows 
us to calculate an H2 column density by subtracting the 
HI column density. This in turn can be compared to the 
CO emission from the Taurus molecular cloud. If the CO to 
H2 conversion factor held down to the smallest scales, then 
the slope of an N-r^ vs. Wco plot would yield this CO X- 
factor. In practice, we would expect this relation to break 
down in resolved clouds because of its local fluctuations. It 
is further known that the X-factor also varies on a cloud-by- 
cloud basis (e.g. Burgh et al. 20071. For these reasons, we 



caution that the comparison between our A^H2 and Wco is 
for the purpose of a self-consistency check only. 

In the left and right panels of Figure [5] we show the 
derived H2 column density vs. Wco for the pixels in the 
map with HINSA strength values of at least 4 and 15, re- 
spectively, using our adopted filament depth of 5 parsec. In 
the left panel we can distinguish roughly three features: a 
main cloud of points, a tail at the bottom (towards the right; 
more about this below) and another tail towards the upper- 
left part of the plot. The latter points disappear entirely for 
the higher HINSA strength threshold. 

We also plot a linear correlation with a slope equal to 
the CO factor derived by Dame et al. (20011 on the ba- 



sis of far-infrared imaging combined with HI maps, as well 
as dashed lines corresponding to twice and half this value, 
respectively. These dashed lines correspond also approxi- 



mately to the reported range of the X-factor (cf. Section 
[2|. We see that the points fall roughly in this range, but a 
least-squares fit of the points themselves, shown by the solid 
gray line, yields an even shallower slope of ~ 0.2 times the 



Dame et al. (20011 X-factor 



To investigate this correlation more closely, we defined 
five regions of interest based on the velocity-integrated local 
CO emission, as shown in Figure |10| These regions trace 
denser parts of the cloud as shown in the right panel of 
Figure [s] Figure 11 (left panel) shows the pixels of these five 
regions in the Nh2-Wco plane with different symbols. For 
each one of the regions, the slopes are shallow and lacking 
any real correlation, but cumulatively a slope is obtained of 
0.9 X 10^° cm"^(K km s~^)~^, which is on the lower end of 
the range of usually quoted X-factors. 

It can be seen that the derived A^H2 for the points in 
region 1, on the eastern tip of the main Taurus area, is con- 
sistently low when compared to the observed CO emission. 
We speculate that our filament depth assumption is off in 
this area. If the filament is curved towards or away from the 
observer, or if the cloud resembles less of a filament in this 
region, then we underestimate the filament depth and as a 
consequence we consistently underestimate the total hydro- 
gen column density. If the cloud is twice as thick here than 
assumed, this would put the points in line with the other 



points. Indeed, Torres et al. (20091 note that Taurus may 



indeed be as deep as it is long in certain places (25 parsecs, 
according to their measurement). 

Region 3 does not have many points left after impos- 
ing a HINSA-strength of 15. Region 4 centers on what is 
thought to be a parallel filament in the Taurus cloud and 
assuming the same filament depth results in column densi- 
ties consistent with those found in the main part. The same 
can be said of region 5, which is at the north-west end. 

Comparing with Figure|9](left panel), we can see that in 
Figure [TT] (left), the points in the upper-left corner have dis- 
appeared - they are not a part of the five selected regions in 
the main Taurus region. One possibility is that these points 
originate from a part of Taurus with a smaller column depth. 
Then our assumed column depth would be too large, lead- 
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Figure 9. Comparison of predicted N}j ^ vs. C O of the full map thresholded by HINSA map values of 4 (left panel) and 15 (right panel). 
The black solid line is the |Dame et al.| ( |2001| l CO X-factor. The dashed lines represent twice and half the solid line's slope. The gray 
solid line shows a simple least squares fit to the data points. The size of the circles is proportional to the HINSA map strength. 




Figure 10. Regions of interest (boxes) displayed on top of a CO 
emission map. 



ing us to overestimate total hydrogen column density and 
therefore also the molecular hydrogen column density. Al- 
ternatively, if the area to which these points belong is part 
of a filamentary structure, and has a column depth consis- 
tent with the one we assumed, the measured A'^hi could be 
too far from the center of that filamentary structure, result- 
ing in an underestimated HI column because of geometric 
considerations, as we showed in Figure [6j and therefore an 
overestimated total hydrogen volume density. 

We show the effect of using a different cloud depth or 
a different x in the right panel of Figure [Tl] a depth of 10 
pc is used here or a x of 10, both yielding the same result 
since n is linearly proportional to both of them. Recalling 
Equation |l| and N = nx ddoud, it follows immediately that 
N oc Godcioud. 

It is important to stress that our results are dependent 
on using an accurate estimate of the incident UV flux Go. 
Since we used one constant value, we would expect that a 



better model of the radiation field would improve the accu- 
racy of our results. Furthermore, we can only calculate an 
H2 column density if we assume a certain cloud depth. 



8 RANDOM MAPS 

Because the distribution of points in Figures [9land[TT] ex- 
hibits significant amounts of scatter, it is useful to com- 
pare those results to cases in which the HI intensities in 
the maps are randomized, in order to see whether our re- 
sults in those figures provide useful information. For this 
purpose, we considered the atomic hydrogen map random- 
ized in two different ways. Then, we proceeded to calculate 
the total hydrogen column densities in the normal fashion, 
to see whether our results contain more information than 
these random cases. 

First, we replaced all the pixels in the field by random 
values in the range of detected HI column densities. Based on 
the Taurus HI map we created a map with random column 
densities ranging from 2.5 x 10^° to 9.3 x 10^" cm . These 
values are uniformly distributed within this range. Secondly, 
we took the observed HI column density map and randomly 
swapped pixel values. This way the global distribution of 
column density values is preserved, but the spatial image is 
destroyed. 

Then, we calculated total hydrogen volume densities 
from both of these maps and compared the calculated H2 
column densities with the real CO emission map (Figure 
121. The point sizes are proportional to the original (non- 
randomized) HINSA strength at each location. It can be 
seen that no clear correlation is left. The left panel plot 
(uniformly distributed Ahi) shows more values towards the 
upper part of the plot than the right panel plot (shuffled 
Ahi) because the distribution of the observed HI map is 
not uniform and contains fewer lower-end Ahi values. For 
comparison, see the same plot for the column densities cal- 
culated using original HI map as was shown in Figure [9] A 
much more restricted range of values is found here as well 
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Figure 11. Comparison of predicted N[]^ vs. CO of the 5 selected regions as marked on Figure [To| A cloud depth of 5 parsec is assumed 
in the left panel. For each region a fitted least-square solid gray line is shown. The dash-dotted gray line is the fitted slope for all regions 
combined. The right panel shows the same plot, but assuming either a cloud depth of 10 parsec or x = 10, since the dependency of A'^ 
on X and cloud depth is the same. The size of the plot symbols is proportional to the HINSA map strength. 




as a better correlation, and therefore we conclude that our 
results provide real information. 



9 COMPARISONS TO FAR-INFRARED 
EMISSION 

Dust emission in the far-infrared has been found to corre- 
late well with the total hydrogen column density. [Dame et al.| 
(20011) used thelSchlegel et al.l (119981) IRAS 100 /.im to calcu- 



late the CO X-factor in our own Galaxy. To this end, they 
calculated the ratio between the atomic hydrogen column 
density and the 100 /im emission in regions free of CO(l-O) 



emission. The resulting ratio, 0.9 x 10 cm 



7(MJy 



(with an estimated uncertainty of 50%) can then be used to 
calculate the total hydrogen column densities in the regions 
with CO(l-O) emission which in turn leads to the calculation 
of the CO X-factor. It should be noted that while iDame et alT] 



("20011 allowed for local variations of this 100 pim conversion 
factor, we adopt their globally averaged ratio here. 



Abergel et al. ( 1995 1 compared the far-infrared emission 



of Taurus to CO emission and found that they correlate well 
globally. They estimate that at least 70% of the 100 ^m- 
traced matter is traced by CO (1-0) emission as well. In other 
words, a bit more gas should be traced by 100 /im emission 
than by CO(l-O) in the Taurus region. For reference, we 
plot the HINSA strength contours on top of the FIR map in 
Figure 13 (cf. Figure |4|. It can easily be seen that the FIR 
map and the CO map are qualitatively very similar in the 
densest regions. However, the infrared map shows additional 
diffuse structure not present in the CO map, consistent with 
the notion that the far-infrared continuum traces the total 
hydrogen gas rather than only the molecular gas. 

We compare our results to the same (temperature- 
corrected) far-infrared map for completeness. A plot of this 
comparison is shown in the left panel of Figure [141 ^ cloud 
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Figure 14. Left panel: Ntot (our method) vs. Ntot (lOO/xm), with a 1-to-l line for comparison using < N[fj/Iioo >= 0.9 X 
10^" cm~^/(MJy sr~i). The gray dots have a HINSA strength between 4 and 15; the black dots have a HINSA strength above 15. 
The latter can be seen to under-predict the total gas column densities as compared to the 100/im-predicted ones. Additionally, the same 
regions as marked in Figure [To] are marked. Right panel: same plot but the Ntot values of regions 1 and 2 have been artificially multiplied 
by 2. 
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Figure 13. Contours of HINSA strength on top of the 
(temperature-corrected) lOO/^m map. HINSA strength contours 
are 4, 12 and 20 K km-^ s^ 



depth of 5 parsecs is assumed as before. The FIR map was 
smoothed and reprojected to the same resolution as the CO 
map, namely 8 arcmin. The gray pixels correspond to a 
HINSA strength of 4 to 15; the black pixels to a HINSA 
strength of 15 and above, and additionally the same regions 
as marked in Figure [lO] are indicated. There is a feature at 
the highest values of the 100 /im emission that is caused by 
a single source in the north-west part of the map (region 
5) that partially overlaps with a HINSA detection. This ex- 
plains the outlying points on the right-hand side of Figure 

M 

It can be seen that the pixels corresponding to the 
higher HINSA strength are mostly confined to a single area 
in the plot (and regions 1 and 2). If we assume that the 



higher HINSA strength corresponds to a better applicability 
of our model, then this suggests that the model is underes- 
timating the total hydrogen column density somewhat (the 
black pixels indicating a HINSA strength of 15 are predomi- 
nantly situated below the 1-to-l correlation line). A possible 
explanation of this deficiency would be that we have under- 
estimated the column depth in region 1 and underestimated 
Go in region 2, both by a factor ~ 2 (note that A'^ is linearly 
proportional to both column depth and Go). The right panel 
of Figure [l4] shows the result of multiplying our by a 
factor of 2 in these regions. A tighter correlation is obtained 
in this case. (See also the discussion of the radiation field in 



Section 10 



Even though the comparison of our method to far- 
infrared emission shows a clear structure, if we conserva- 
tively estimate our results to have an uncertainty compara- 
ble to the infrared method, namely 50%, the correlation is 
significantly degraded. We expect the derived total hydro- 
gen volume densities to have the same or better uncertainties 
than the ones in the nearest extragalactic case, namely M33 
( [Heiner et aL][20TT| ). The same range of values is obtained 
with our method as with the far-infrared method, meaning 
that both methods yield consistent results. In order to ob- 
tain better statistics a wider range of molecular clouds will 
need to be investigated in a similar fashion, both lower den- 
sity as well as higher density clouds up to the applicability 
limit of ~ 10* cm~^ total hydrogen density. 



10 CHARACTERIZING THE RADIATION 
FIELD IN TAURUS 

Considering the relatively low implied H2 column densities 
that we found in Taurus, we will discuss in this section the 
possibility of a stronger UV radiation field than we assumed. 
A stronger field would increase our inferred H2 column den- 
sities. 
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10.1 PDR model comparison 



Allen et al. (20041 presented calculations from a model of 



a one-dimensional plane-parallel PDR and how CO and HI 
emission relate to the incident flux Go and the total hydro- 
gen volume density: Their Figures 1 and 2 provide diagnostic 
diagrams predicting a limited range of Go and n given an 
HI column density and CO intensity. Their work was an ex- 
tension of the modeling done by Kaufman et al.| ( |1999 1 that 
solved the chemistry, radiative transfer and thermal balance 
of one-dimensional PDRs. The , Allen et al.| (|2004) results 
predict that the CO emission at the expected density of the 
Taurus region is in the regime where this emission is only 
very weakly dependent on the incident UV flux. [Alien et al.| 
( 2004 1 compared their calculations to the overall character- 



istics of G216-2.5 (Maddalena's Cloud). We will now do the 
same for the Taurus area. 

Taking the average CO emission over the full map yields 
(Wco) = 20.3 K km s'^ The average HI column density in 
the Taurus area in our selected velocity range is (A^hi) = 
K 10^0 cm"^. Intersecting the values of {Wco) and (Ahi) 
Allen et al. (20041, their Figure 2, implies Go ~ 6 and 



a total hydrogen volume density of several hundred cm~^. 
Computing the HI column density over a broader velocity 
range would raise the HI column density, implying a higher 
incident UV flux if n is fixed. 



Figure 20 in HoUenbach et al. ( 1991 1 provides an ad- 



ditional way to infer Go from the total hydrogen volume 
density and the ratio Wco/^ioo- Taking the average 100 ^.m 
emission in the Taurus area map of 40 MJy sr~^ we obtain 
{Wco) / {I wo) = 0.5, which for volume densities of several 
hundred cm~^ implies Go > 10. 

The above discussion suggests that Go is in the range 
of 1 to perhaps several times 10. We consider this quite 
realistic, considering the extragalactic case - for example 



see the range of Go values obtained in M33 (Heiner et al 



2011). All that is needed is a cluster of a few dozen O or 



B stars within a few hundred parsecs distance (further than 
our own distance from Taurus) and a sufficiently clumpy 
interstellar medium to let the UV radiation penetrate. 



10.2 Far-ultraviolet emission 

A far-ultraviolet map of the Taurus region was presented 



by Lee et al. (20061, although the emission is comparatively 



faint (x 0.2). They argue that the observed FUV intensity 
comes mainly from foreground radiation while radiation at 
the back of the cloud is blocked. Far-infrared emission should 
not be blocked, potentially providing an alternative measure 
of the UV field, but they point out that the relationship 
between /puv and Iioo/jm should deviate from a simple linear 
relationship in optically thick regions (such as Taurus). Such 
a linear relationship was proposed by Haikala et al. (19951, 



who provide their own estimate of the interstellar radiation 
field (see also |Murthy fc Sahnow]|2004 l. 

It is still instructive to convert the 100 /.tm image to 
an equivalent map of x (the far-ultraviolet fiux scaled by 
their interstellar radiation field value) using the linear rela- 
tion from Haikala et al. ( 1995 1. The conversion with scaling 
becomes: 



where x is in equivalent dimensionless units of the usual x 
since it was scaled to the interstellar radiation field value in 
the solar neighborhood. Scaling the 100 /im image in this 
fashion yields values up to x = 3, below our assumed x = 
5. However, Figure 6 of Murthy & Sahnow (20041 shows 



the correlation between 7fuv and /loo for a broader range 
of values. It can be see that the inferred values of x can 
be several times higher because of significant scatter in the 
correlation for values similar to the ones found in the Taurus 



11 CONCLUSIONS 

We used a simple one-dimensional PDR model to predict to- 
tal hydrogen volume densities from HI column densities in 
the Taurus molecular cloud in regions with strong HI nar- 



row self-absorption (HINSA). As Allen et al. (20041 noted 



this procedure is primarily of use in the extragalactic case. 
However, it is important to test the method by applying it 
locally in our own Galaxy in order to validate the extra- 
galactic results. To this end, we selected atomic hydrogen 
column densities from a limited range of velocities around 
the main molecular cloud radial velocity (ignoring HINSA 
corrections to the column densities). 

The obtained total hydrogen volume densities appear to 
follow a log-normal distribution, potentially with a power- 
law tail at the high density end. Our method provides a 
promising tool to obtain the probability density function 
of the volume density of the gas in molecular clouds. We 
find an average total hydrogen volume density of 436 cm""^ 
(cr = 129 cm~^) on the basis of pixels selected for HINSA 
features. 

For the first time, we make a direct detailed com- 
parison between our PDR-based method and CO-derived 
column densities (cf. [Allen et al.| [2004 [ ) as well as far- 
infrared-derived column densities. We found a correlation 
between A''h2 and Wco with a cumulative slope of 0.9 x 
10^° cm~^(K km s~^)~^, implying a relatively low value of 
the CO to H2 conversion factor, at least when the full en- 
semble of data points is considered (individual parts of the 
Taurus cloud essentially show no correlation. However, as- 
suming that the far-infrared emission is a good tracer of 
the total gas, our method yields lower total hydrogen col- 
umn densities as compared to the column densities implied 
by the far-infrared emission, at least for higher values of the 
HINSA strength. There are three ways in which our total hy- 
drogen column densities can be made larger: using a larger 
cloud depth, a higher incident UV flux Go, or a lower HI 
column density (obtained by narrowing the velocity range 
of atomic hydrogen associated to Taurus). 



We adopted x = 5 (Co = 4.25) from Goldsmith et al 



(20101, although in line with their findings, this seems to 



X = (128 /loo - 264) / 6490, 



(2) 



be the lowest value for the UV radiation field that fits the 
observations. Twice this value might yield better results. We 
conclude, based on our results and the available literature, 
that a higher Go of at least 10 in the general Taurus area 
would be expected. 

In the case of the central region of Taurus, it appears 
that all the gas mass is accounted for (since our results 
are consistent with both CO and FIR observations), unless 
the true Go pervading the molecular cloud is, unexpectedly, 
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at least several times higher than previously assumed. We 
might expect to find hidden 'dark' gas based on the work 
by Wolfire et al. (20101, or the observations of galactic su- 
pershells by Dawson et al. (20111. The parts of the Tau- 



rus area where we applied the plane-parallel PDR model 
were not CO-selected, but rather by the presence of HINSA; 
however, all selected areas do have detectable CO emission. 
One could assume that molecular gas not traced by CO is 
present in the atomic hydrogen layer surrounding the CO 
cloud, but the results of our method show no evidence of 
this in the form of predicted total hydrogen column densi- 
ties consistently higher than those inferred by CO emission 
or infrared emission. Molecular gas not traced by CO could 
still be present in the areas directly surrounding Taurus, 
further away from the CO cloud on the plane of the sky. 
However, the lack of better information about the detailed 
UV radiation field in the Taurus region prevents a thorough 
search at this time. It would be of interest to apply the same 
model to similar clouds in our Galaxy with hints of 'dark' 
gas that is not traced by CO or where higher X-factors have 
been inferred. 

Assuming a consistent relation between our inferred 
molecular gas versus the CO-implied gas, we speculate that 
the eastern tip of the Taurus cloud has twice the thickness 
of the main area, independent of what thickness is adopted 
or will be found for the main area. 

The results of this work could be improved through a 
combination of cloud depth information to aid the direct 
comparison to CO (anticipating e.g. the Gould Belt Dis- 
tance Survey, Loinard et al.|20"TT I and a more sophisticated 
way to calculate local (related) HI column densities. More 
accurate HI column density measurements would require, for 
example, a way to distinguish foreground and background 
HI emission from emission from the molecular cloud, if we 
could find out more about the local velocity and tempera- 
ture structure. 

Going forward, we expect to be able to use the simple 
PDR model in future work as a tool for investigating the 
atomic-molecular transition in molecular clouds. 
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